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Abstract 

Inclusive jet differential cross sections have been measured in neutral current 
deep inelastic e+p scattering for boson virtualities > 125 GeV^. The data 
were taken using the ZEUS detector at HERA and correspond to an integrated 
luminosity of 38.6 pb~^. Jets were identified in the Breit frame using the longi- 
tudinally invariant kx cluster algorithm. Measurements of differential inclusive 
jet cross sections are presented as functions of jet transverse energy (E'^j^^J, jet 
pseudorapidity and Q^, for jets with -Exjct > 8 GeV. Next-to-leading-order QCD 
calculations agree well with the measurements both at high and high E^-^^. 
The value of as{Mz), determined from an analysis of da/dQ^ for > 500 GeV^, 
is as{Mz) = 0.1212 ±0.0017 (stat.)+J:S (syst.)^!;:S (th.). 



The ZEUS Collaboration 

S. Chekanov, D. Krakauer, S. Magill, B. Musgrave, J. Repond, R. Yoshida 
Argonne National Laboratory, Argonne, Illinois 60439-4815 

M.C.K. Mattingly 

Andrews University, Berrien Springs, Michigan 49104-0380 

P. Antonioli, G. Bari, M. Basile, L. Bellagamba, D. Boscherini, A. Bruni, G. Bruni, 
G. Cara Romeo, L. Cifarelli, F. Cindolo, A. Contin, M. Corradi, S. De Pasquale, P. Giusti, 
G. lacobucci, A. Margotti, R. Nania, F. Palmonari, A. Pesci, G. Sartorelli, A. Zichichi 
University and INFN Bologna, Bologna, Italy ^ 

G. Aghuzumtsyan, D. Bartsch, I. Brock, J. Crittenden^, S. Goers, H. Hartmann, E. Hilger, 

P. Irrgang, H.-P. Jakob, A. Kappes, U.F. Katz^, R. Kerger^, O. Kind, E. Paul, J. Rautenberg^, 
R. Renner, H. Schnurbusch, A. Stifutkin, J. Tandler, K.C. Voss, A. Weber 
Physikalisches Institut der Universitdt Bonn, Bonn, Germany * 

D. S. Bailey^ N.H. Brook^, J.E. Cole, B. Foster, G.P. Heath, H.F. Heath, S. Robins, 

E. Rodrigues^ J. Scott, R.J. Tapper, M. Wing 

H. H. Wills Physics Laboratory, University of Bristol, Bristol, United Kingdom^ 

M. Capua, A. Mastroberardino, M. Schioppa, G. Susinno 

Calabria University, Physics Department and INFN, Cosenza, Italy 

J.Y. Kim, Y.K. Kim, J.H. Lee, IT. Lim, M.Y. Pac^ 
Chonnam National University, Kwangju, Korea ^ 

A. Caldwell^ M. Helbich, X. Liu, B. Mellado, Y. Ning, S. Paganis, Z. Ren, W.B. Schmidke, 

F. Sciulh 

Nevis Laboratories, Columbia University, Irvington on Hudson, New York 10027° 

J. Chwastowski, A. Eskreys, J. Figiel, K. Olkiewicz, K. Piotrzkowski^, M.B. Przybycieh^°, 
P. Stopa, L. Zawiejski 

Institute of Nuclear Physics, Cracow, Poland * 

L. Adamczyk, T. Bold, I. Grabowska-Bold, D. Kisielewska, A.M. Kowal, M. Kowal, 
T. Kowalski, M. Przybycieh, L. Suszycki, D. Szuba, J. Szuba^^ 

Faculty of Physics and Nuclear Techniques, University of Mining and Metallurgy, Cracow, 
Poland ^ 

A. Kotahski^^, W. Slomihski^^ 

Department of Physics, Jagellonian University, Cracow, Poland 



I 



L.A.T. Bauerdick^^ , U. Behrens, K. Borras, V. Chiochia, D. Dannheim, M. Derrick^^, 
G. Drews, J. Fourletova, A. Fox-Murphy, U. Fricke, A. Geiser, F. Goebel*, P. Gottlicher^^, 

0. Gutsche, T. Haas, W. Hain, G.F. Hartner, S. Hillert, U. Kotz, H. Kowalski^^, G. Kram- 
berger, H. Labes, D. Lelas, B. Lohr, R. Mankel, M. Martmez^^, I. -A. Melzer-Pellmann, 
M. Moritz, D. Notz, M.C. Petrucci^^ A. Polini, A. Raval, U. Schneekloth, F. Selonke^^ 
B. Surrow^'', H. Wessoleck, R. Wichmann^^, G. Wolf, C. Youngman, W. Zeuner 
Deutsches Elektronen- Synchrotron DESY, Hamburg, Germany 

A. Lopez-Duran Viani^^, A. Meyer, S. Schlenstedt 
DESY Zeuthen, Zeuthen, Germany 

G. Barbagli, E. Gallo, C. Genta, P. G. Pelfer 
University and IN FN, Florence, Italy ^ 

A. Bamberger, A. Benen, N. Coppola, H. Raach 

Fakultdt fiir Physik der Universitdt Freiburg i.Br., Freiburg i.Br., Germany * 

M. Bell, P.J. Bussey, A.T. Doyle, C. Glasman, S. Hanlon, S.W. Lee, A. Lupi, G.J. Mc- 
Cance, D.H. Saxon, I.O. Skillicorn 

Department of Physics and Astronomy, University of Glasgow, Glasgow, United King- 
dom 

1. Gialas 

Department of Engineering in Management and Finance, Univ. of Aegean, Greece 

B. Bodmann, T. Carli, U. Holm, K. Klimek, N. Krumnack, E. Lohrmann, M. Milite, 

H. Salehi, S. Stonjek^^, K. Wick, A. Ziegler, Ar. Ziegler 

Hamburg University, Institute of Exp. Physics, Hamburg, Germany ^ 

C. Collins-Tooth, C. Fondas, R. Gongalo^ K.R. Long, F. Metlica, D.B. Miller, A.D. Tap- 
per, R. Walker 

Imperial Gollege London, High Energy Nuclear Physics Group, London, United King- 
dom ™ 

P. Cloth, D. Filges 

Forschungszentrum Jiilich, Institut fiir Kernphysik, Jiilich, Germany 

M. Kuze, K. Nagano, K. Tokushuku^'' , S. Yamada, Y. Yamazaki 
Institute of Particle and Nuclear Studies, KEK, Tsukuba, Japan ^ 

A.N. Barakbaev, E.G. Boos, N.S. Pokrovskiy, B.O. Zhautykov 

Institute of Physics and Technology of Ministry of Education and Science of Kazakhstan, 
Almaty, Kazakhstan 

H. Lim, D. Son 

Kyungpook National University, Taegu, Korea ^ 



II 



F. Barreiro, O. Gonzalez, L. Labarga, J. del Peso, I. Redondo^^, J. Terron, M. Vazquez 
Departamento de Fisica Teorica, Universidad Autonoma de Madrid, Madrid, Spain ' 

M. Barbi, A. Bertolin, F. Corriveau, A. Ochs, S. Padhi, D.G. Stairs, M. St-Laurent 
Department of Physics, McGill University, Montreal, Quebec, Canada H3 A 2T8 

T. Tsurugai 

Meiji Gakuin University, Faculty of General Education, Yokohama, Japan 

A. Antonov, P. Danilov, B.A. Dolgoshein, D. Gladkov, V. Sosnovtsev, S. Suchkov 
Moscow Engineering Physics Institute, Moscow, Russia ^ 

R.K. Dementiev, P.F. Ermolov, Yu.A. Golubkov, I.I. Katkov, L.A. Khein, I. A. Korzhav- 
ina, V.A. Kuzmin, B.B. Levchenko, O.Yu. Lukina, A.S. Proskuryakov, L.M. Shcheglova, 
N.N. Vlasov, S.A. Zotkin 

Moscow State University, Institute of Nuclear Physics, Moscow, Russia ^ 

C. Bokel, J. Engelen, S. Grijpink, E. Koffeman, P. Kooijman, E. Maddox, A. Pellegrino, 
S. Schagen, E. Tassi, H. Tiecke, N. Tuning, J.J. Velthuis, L. Wiggers, E. de Wolf 
NIKHEF and University of Amsterdam, Amsterdam, Netherlands ^ 

N. Briimmer, B. Bylsma, L.S. Durkin, J. Gilmore, CM. Ginsburg, C.L. Kim, T.Y. Ling 

Physics Department, Ohio State University, Golumbus, Ohio 43210"" 

S. Boogert, A.M. Cooper-Sarkar, R.C.E. Devenish, J. Ferrando, G. Grzelak, T. Mat- 
sushita, M. Rigby, O. Ruske^^, M.R. Sutton, R. Walczak 
Department of Physics, University of Oxford, Oxford United Kingdom 

R. Brugnera, R. Carlin, F. Dal Corso, S. Dusini, A. Garfagnini, S. Limentani, A. Longhin, 

A. Parenti, M. Posocco, L. Stanco, M. Turcato 

Dipartimento di Fisica dell' Universitd and INFN, Padova, Italy ^ 

E.A. Heaphy, B.Y. Oh, P.R.B. SaulP^ J.J. Whitmore^^ 

Department of Physics, Pennsylvania State University, University Park, Pennsylvania 
16802° 

Y. Iga 

Polytechnic University, Sagamihara, Japan ^ 

G. D'Agostini, G. Marini, A. Nigro 

Dipartimento di Fisica, Universitd 'La Sapienza' and INFN, Rome, Italy ^ 
C. Cormack29, J.C. Hart, N.A. McCubbin 

Rutherford Appleton Laboratory, Ghilton, Didcot, Oxon, United Kingdom ™ 
C. Heusch 

University of Galifornia, Santa Gruz, Galifornia 95064 " 



III 



I.H. Park 

Department of Physics, Ewha Womans University, Seoul, Korea 
N. Pavel 

Fachbereich Physik der Universitdt- Gesamthochschule Siegen, Germany 
H. Abramowicz, A. Gabareen, S. Kananov, A. Kreisel, A. Levy 

Raymond and Beverly Sackler Faculty of Exact Sciences, School of Physics, Tel-Aviv 
University, Tel- Aviv, Israel 

T. Abe, T. Fusayasu, S. Kagawa, T. Kohno, T. Tawara, T. Yamashita 
Department of Physics, University of Tokyo, Tokyo, Japan ^ 

R. Hamatsu, T. Hirose^^, M. Inuzuka, S. Kitamura'^'^, K. Matsuzawa, T. Nishimura 
Tokyo Metropolitan University, Deptartment of Physics, Tokyo, Japan ^ 

M. Arneodo^^ N. Cartiglia, R. Cirio, M. Costa, M.I. Ferrero, S. Maselli, V. Monaco, 

C. Peroni, M. Ruspa, R. Sacchi, A. Solano, A. Staiano 

Universitd di Torino, Dipartimento di Fisica Sperimentale and INFN, Torino, Italy ^ 

R. Galea, T. Koop, G.M. Levman, J.F. Martin, A. Mirea, A. Sabetfakhri 
Department of Physics, University of Toronto, Toronto, Ontario, Canada M5S 1A7 

J.M. Butterworth, C. Gwenlan, R. Hall-Wilton, T.W. Jones, M.S. Lightwood, J.H. Loizides^^^ 
B.J. West 

Physics and Astronomy Department, University College London, London, United King- 
dom 

J. Ciborowski^s, R. Ciesielski^^ , R.J. Nowak, J.M. Pawlak, B. Smalska3^ J. Sztuk3^ 

T. Tymieniecka^^ , A. Ukleja^^, J. Ukleja, A.F. Zarnecki 

Warsaw University, Institute of Experimental Physics, Warsaw, Poland ^ 

M. Adamus, P. Plucinski 

Institute for Nuclear Studies, Warsaw, Poland 

Y. Eisenberg, L.K. Gladilin'^^, D. Hochman, U. Karshon 

Department of Particle Physics, Weizmann Institute, Rehovot, Israel ^ 

D. Kgira, S. Lammers, L. Li, D.D. Reeder, A. A. Savin, W.H. Smith 
Department of Physics, University of Wisconsin, Madison, Wisconsin 53706 

A. Deshpande, S. Dhawan, V.W. Hughes, P.B. Straub 

Department of Physics, Yale University, New Haven, Connecticut 06520-8121 " 

S. Bhadra, CD. Catterall, S. Fourletov, S. Menary, M. Soares, J. Standage 
Department of Physics, York University, Ontario, Canada MS J IPS " 



IV 



^ now at Cornell University, Ithaca/NY, USA 

^ on leave of absence at University of Erlangen-Niirnberg, Germany 
^ now at Ministere de la Culture, de L'Enseignement Superieur et de la Recherche, Lux- 
embourg 

^ supported by the GIF, contract 1-523-13.7/97 
5 PPARC Advanced fellow 

^ supported by the Portuguese Foundation for Science and Technology (FCT) 
^ now at Dongshin University, Naju, Korea 
* now at Max-Planck- Institut fiir Physik, Miinchen/ Germany 
^ now at Universite Catholique de Louvain, Louvain-la-Neuve/Belgium 
now at Northwestern Univ., Evanston/IL, USA 

partly supported by the Israel Science Foundation and the Israel Ministry of Science 
supported by the Polish State Committee for Scientific Research, grant no. 2 P03B 
09322 

member of Dept. of Computer Science 
now at Fermilab, Batavia/IL, USA 
on leave from Argonne National Laboratory, USA 
16 now at DESY group FEB 

1^ on leave of absence at Columbia Univ., Nevis Labs., N.Y./USA 
1^ now at INFN Perugia, Perugia, Italy 
1^ retired 

"^^ now at Brookhaven National Lab., Upton/NY, USA 
^1 now at Mobil com AG, Rendsburg-Biidelsdorf, Germany 

now at Deutsche Borse Systems AG, Frankfurt/Main, Germany 
^■^ now at Univ. of Oxford, Oxford/UK 

also at University of Tokyo 

now at LPNHE Ecole Polytechnique, Paris, France 
^6 now at IBM Global Services, Frankfurt/Main, Germany 
^'^ now at National Research Council, Ottawa/Canada 

on leave of absence at The National Science Foundation, Arlington, VA/USA 

now at Univ. of London, Queen Mary College, London, UK 

present address: Tokyo Metropolitan University of Health Sciences, Tokyo 116-8551, 
Japan 

^1 also at Universita del Piemonte Orientale, Novara, Italy 
supported by Argonne National Laboratory, USA 
also at Lodz University, Poland 

supported by the Polish State Committee for Scientific Research, grant no. 2 P03B 
07222 

now at The Boston Consulting Group, Warsaw, Poland 



V 



Lodz University, Poland 

supported by German Federal Ministry for Education and Research (BMBF), POL 
01/043 

on leave from MSU, partly supported by University of Wisconsin via the U.S. -Israel BSF 



VI 



supported by the Natural Sciences and Engineering Research Council of 
Canada (NSERC) 

supported by the German Federal Ministry for Education and Research 
(BMBF), under contract numbers HZIGUA 2, HZIGUB 0, HZIPDA 5, 
HZIVFA 5 

supported by the MINERVA Gesellschaft fiir Forschung GmbH, the Israel Sci- 
ence Foundation, the U.S. -Israel Binational Science Foundation, the Israel Min- 
istry of Science and the Benozyio Center for High Energy Physics 
supported by the German-Israeli Foundation, the Israel Science Foundation, 
and by the Israel Ministry of Science 

supported by the Italian National Institute for Nuclear Physics (INFN) 
supported by the Japanese Ministry of Education, Science and Culture (the 
Monbusho) and its grants for Scientific Research 

supported by the Korean Ministry of Education and Korea Science and Engi- 
neering Foundation 

supported by the Netherlands Foundation for Research on Matter (FOM) 
supported by the Polish State Committee for Scientific Research, grant no. 
620/E-77/SPUB-M/DESY/P-03/DZ 247/2000-2002 

partially supported by the German Federal Ministry for Education and Re- 
search (BMBF) 

supported by the Fund for Fundamental Research of Russian Ministry for 
Science and Education and by the German Federal Ministry for Education 
and Research (BMBF) 

supported by the Spanish Ministry of Education and Science through funds 
provided by CICYT 

supported by the Particle Physics and Astronomy Research Council, UK 
supported by the US Department of Energy 
supported by the US National Science Foundation 

supported by the Polish State Committee for Scientific Research, grant no. 
112/E-356/SPUB-M/DESY/P-03/DZ 301/2000-2002, 2 P03B 13922 
supported by the Polish State Committee for Scientific Research, grant no. 
115/E-343/SPUB-M/DESY/P-03/DZ 121/2001-2002, 2 P03B 07022 



VII 



1 Introduction 



Jet production in neutral current deep inelastic e~^p scattering at high (where is the 
negative of the square of the virtuality of the exchanged boson) provides a testing ground 
for the theory of the strong interaction between quarks and gluons, namely quantum 
chromodynamics (QCD). In deep inelastic scattering (DIS), the predictions of perturbative 
QCD (pQCD) have the form of a convolution of matrix elements with parton distribution 
functions (PDFs) of the target hadron. The matrix elements describe the short-distance 
structure of the interaction and are calculable in pQCD at each order, whilst the PDFs 
contain the description of the long-distance structure of the target hadron. 

The evolution of the PDFs with the scale at which they are probed is predicted in pQCD 
to follow a set of renormalisation group equations (DGLAP equations [1]). However, an 
explicit determination of the PDFs requires experimental input. A wealth of data from 
fixed-target [2] and collider [3,4] experiments has allowed an accurate determination of the 
proton PDFs [5-10]. Good knowledge of PDFs makes measurements of jet production 
in DIS a sensitive test of the pQCD predictions of the short- distance structure of the 
partonic interactions. 

The hadronic final state in neutral current DIS may consist of jets of high transverse 
energy produced in the short-distance process as well as the remnant (beam jet) of the 
incoming proton. A jet algorithm should distinguish as clearly as possible between the 
beam jet and the hard jets. Working in the Breit frame [11] is preferred, since it provides 
a maximal separation between the products of the beam fragmentation and the hard jets. 
In this frame, the exchanged virtual boson (V^*, with V* = 7, Z) is purely space-like, with 
3-momentum q = (0,0, —Q). In the Born process, the virtual boson is absorbed by the 
struck quark, which is back-scattered with zero transverse momentum with respect to the 
V* direction, whereas the beam jet follows the direction of the initial struck quark. Thus, 
the contribution due to the current jet in events from the Born process is suppressed by 
requiring the production of jets with high transverse energy in this frame. Jet production 
in the Breit frame is, therefore, directly sensitive to hard QCD processes, thus allowing 
direct tests of the pQCD predictions. The use of the cluster algorithm [12] to define 
jets in the Breit frame facilitates the separation of the beam fragmentation and the hard 
process in the calculations [13]. 

At leading order (LO) in the strong coupling constant, a^, the boson-gluon-fusion (BGF, 
V*g qq) and QCD-Compton (QCDC, V*q — > qg) processes give rise to two hard jets 
with opposite transverse momenta. The calculation of dijet cross sections in pQCD at 
fixed order in is hampered by infrared-sensitive regions, so that additional jet-selection 
criteria must be applied to make reliable predictions [14]. This complication is absent in 
the case of cross-section calculations for inclusive jet production. 
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This paper presents measurements of several differential cross sections for the inclusive 
production of jets with high transverse energy in the Breit frame. The analysis is re- 
stricted to large values of Q^, > 125 GeV^, and the jets were selected according to 
their transverse energies and pseudorapidities in the Breit frame; in the definition of the 
cross sections, no cut was applied to the jets in the laboratory frame. The measurements 
are compared to next-to-leading-order (NLO) QCD calculations [15] using currently avail- 
able parameterisations of the proton PDFs. The jet selection used allows a reduction in 
the theoretical uncertainty of the NLO QCD calculations with respect to those of dijet 
production [16,17]. A QCD analysis of the inclusive jet cross sections has been per- 
formed, which yields a more precise determination of than was previously possible at 
HERA [17-21]. 



2 Experimental set-up 

The data sample used in this analysis was collected with the ZEUS detector at HERA 
and corresponds to an integrated luminosity of 38.6 ±0.6 pb~^. During 1996-1997, HERA 
operated with protons of energy Ep = 820 GeV and positrons of energy Ee = 27.5 GeV. 
The ZEUS detector is described in detail elsewhere [22,23]. The main components used 
in the present analysis are the central tracking detector [24], positioned in a 1.43 T 
solenoidal magnetic field, and the uranium-scintillator sampling calorimeter (CAL) [25]. 
The tracking detector was used to establish an interaction vertex. The CAL covers 99.7% 
of the total solid angle. It is divided into three parts with a corresponding division 
in the polar angleQ, 6, as viewed from the nominal interaction point: forward (FCAL, 
2.6° < e < 36.7°), barrel (BCAL, 36.7° < 9 < 129.1°), and rear (RCAL, 129.1° < 9 < 
176.2°). The smallest subdivision of the CAL is called a cell. Under test-beam conditions, 
the CAL relative energy resolution is 18%/^E{GeV) for electrons and 35%/^E{GeV) 
for hadrons. Jet energies were corrected for the energy lost in inactive material, typically 
about 1 radiation length, in front of the CAL. The effects of uranium noise were minimised 
by discarding cells in the inner (electromagnetic) or outer (hadronic) sections if they had 
energy deposits of less than 60 MeV or 110 MeV, respectively. A three-level trigger was 
used to select events online [23]. 

The luminosity was measured using the Bethe-Heitler reaction e~^p — > e'^'jp [26]. The 
resulting small-angle energetic photons were measured by the luminosity monitor, a lead- 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction", and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity is 
defined as 77 = — ln(tan |). 
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scintillator calorimeter placed in the HERA tunnel a.t Z = 



-107 m. 



3 Data selection and jet search 

Neutral current DIS events were selected offline using criteria similar to those reported 
previously [27]. The main steps are briefly discussed below. 

The scattered-positron candidate was identifled from the pattern of energy deposits in the 
CAL [28]. The energy [E'^) and polar angle (^e) of the positron candidate were determined 
from the CAL measurements. The variable was reconstructed from the double angle 
method (Q^a) [29], which uses 9e and an angle 7 that corresponds, in the quark-parton 
model, to the direction of the scattered quark. The angle 7 was reconstructed from the 
CAL measurements of the hadronic flnal state [29]. The following requirements were 
imposed on the data sample: 

• a positron candidate of energy E'^ > 10 GeV. This cut ensured a high and well un- 
derstood positron-flnding efflciency and suppressed background from photoproduction 
events, in which the scattered positron escapes down the rear beampipe; 

• He < 0.95, where i/e = I — E'^{1 — cos9e)/{2Ee). This condition removed events in 
which fake positron candidates were found in the FCAL; 

• the total energy not associated with the positron candidate within a cone of radius 0.7 
units in the pseudorapidity-azimuth (77 — 0) plane around the positron direction should 
be less than 10% of the positron energy. This condition removed photoproduction and 
DIS events in which part of a jet was falsely identifled as the scattered positron; 

• for 30° < 6e < 140°, the fraction of the positron energy within a cone of radius 0.3 
units in the r] — (f) plane around the positron direction should be larger than 0.9; for 
6e < 30°, the cut was raised to 0.98. This condition removed events in which a jet was 
falsely identifled as the scattered positron; 

• the vertex position along the beam axis should be in the range —38 < Z < 32 cm; 

• 38 < {E — pz) < 65 GeV, where E is the total energy as measured by the CAL, 
E = Ei, and pz is the Z-component of the vector p = -Ej^i; in both cases the 
sum runs over all CAL cells, Ei is the energy of the CAL cell i and rj is a unit vector 
along the line joining the reconstructed vertex and the geometric centre of the cell 
i. This cut removed events with large initial-state radiation and further reduced the 
background from photoproduction; 

• ^t/VE^ < 2.5 GeV^/^ where is the missing transverse momentum as measured 
with the CAL {"^t = \/v\^ Vy) aiid Et is the total transverse energy in the CAL. 
This cut removed cosmic rays and beam-related background; 



3 



• no second positron candidate with energy above 10 GeV and energy in the GAL, 
after subtracting that of the two positron candidates, below 4 GeV. This requirement 
removed elastic Gompton scattering events (ep ^ip)', 

• Ql^ > 125 GeV2; 

• —0.7 < cos 7 < 0.5. The lower limit avoided a region with limited acceptance due to 
the requirement on the energy of the scattered positron, whilst the upper limit was 
chosen to ensure good reconstruction of the jets in the Breit frame. 

The longitudinally invariant cluster algorithm [12] was used in the inclusive mode [30] 
to reconstruct jets in the hadronic final state both in data and in Monte Garlo (MG) 
simulated events (see Section H). In data, the algorithm was applied to the energy deposits 
measured in the GAL cells after excluding those associated with the scattered-positron 
candidate. The jet search was performed in the pseudorapidity (?7'^)-azimuth {(p^) plane 
of the Breit frame. In the following discussion, E^^^ denotes the transverse energy, rjf the 
pseudorapidity and 0f the azimuthal angle of object i in the Breit frame. For each pair of 
objects (where the initial objects are the energy deposits in the GAL cells), the quantity 

d^J = livf - riff + (0f - <Pff] ■ min(E|?,„ Ef^f (1) 

was calculated. For each individual object, the quantity di = {E^^f was also calculated. 
If, of all the values {dij, d^}, dki was the smallest, then objects k and / were combined into a 
single new object. If, however, dk was the smallest, then object k was considered a jet and 
was removed from the sample. The procedure was repeated until all objects were assigned 
to jets. The jet variables were defined according to the Snowmass convention [31]: 

_ B _ ^T,iVf _ ^T,i^f /^N 

This prescription was also used to determine the variables of the intermediate objects. 

After reconstructing the jet variables in the Breit frame, the massless four-momenta were 
boosted into the laboratory frame, where the transverse energy (E'^^j^J, the pseudorapid- 
ity (?7j^t) and the azimuthal angle (0j^t) of ^ach jet were calculated. Energy corrections 
were then applied to the jets in the laboratory frame and propagated into the jet trans- 
verse energies in the Breit frame. In addition, the jet variables in the laboratory frame 
were used to apply additional cuts on the selected sample: 

• events were removed from the sample if the distance of any of the jets to the positron 
candidate in the rj — (p plane of the laboratory frame. 



was smaller than 1 unit. This requirement removed some background from photopro- 
duction and improved the purity of the sample; 
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• events were removed from the sample if any of the jets was in the backward region of 
the detector (r/j^^. < —2). This requirement removed events in which a radiated photon 
from the positron was misidentified as a hadronic jet in the Breit frame; 

• jets with low transverse energy in the laboratory frame {E^^^^ < 2.5 GeV) were not 
included in the final sample; this cut removed a small number of jets for which the 
uncertainty on the energy correction was large. 

It should be noted that these cuts were applied to improve the efficiency and purity 
of the sample of jets and were not used to define the phase-space region of the cross- 
section measurements. The simulated events were used to correct these effects on the 
cross sections. In particular, the effects of the last two cuts were estimated to be smaller 
than 3%. The final data sample contained 8523 events with at least one jet satisfying 
E^j^t > 8 GeV and -2 < r]f^^ < 1.8. With the above criteria, 5073 one-jet, 3262 two-jet, 
182 three-jet and 6 four-jet events were found. Since the net transverse momentum of 
the hadronic final state in the Breit frame is zero, an event with a single jet, according 
to a given selection criterion, must contain at least one other jet balancing its transverse 
momentum; however, this jet will not necessarily satisfy the jet-selection criteria. 



4 Monte Carlo simulation 

Samples of events were generated to determine the response of the detector to jets of 
hadrons and the correction factors necessary to obtain the hadron-level jet cross sections. 
The generated events were passed through the GEANT 3.13-based [32] ZEUS detector- 
and trigger-simulation programs [23] . They were reconstructed and analysed by the same 
program chain as the data. 

Neutral current DIS events were generated using the LEPTO 6.5 program [33] interfaced 
to HERACLES 4.5.2 [34] via DJANGO 6.2.4 [35]. The HERACLES program includes 
photon and Z exchanges and first-order electroweak radiative corrections. The QCD 
cascade was modelled with the colour-dipole model [36] by using the ARIADNE 4.08 
program [37] and including the BGF process. The colour-dipole model treats gluons 
emitted from quark-antiquark (diquark) pairs as radiation from a colour dipole between 
two partons. This results in partons that are not ordered in their transverse momenta. 
The CTEQ4D [5] proton PDFs were used. As an alternative, samples of events were 
generated using the model of LEPTO based on first-order QCD matrix elements plus 
parton showers (MEPS). For the generation of the samples with MEPS, the option for 
soft-colour interactions was switched off [38]. In both cases, fragmentation into hadrons 
was performed using the LUND [39] string model as implemented in JETSET 7.4 [40]. 
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The jet search was performed on the MC events using the energy measured in the CAL 
cells in the same way as for the data. Using the sample of events generated with either 
ARIADNE or LEPTO-MEPS and after applying the same offline selection as for the data, 
a good description of the measured distributions for the kinematic and jet variables was 
found. The same jet algorithm was also applied to the hadrons (partons) to obtain the 
predictions at the hadron (parton) level. The MC programs were used to correct the 
measured cross sections for QED radiative effects. 

5 NLO QCD calculations 

The measurements were compared with NLO QCD {0{a^)) calculations obtained using 
the program DISENT [15]. The calculations were performed in the MS renormalisation 
and factorisation schemes using a generalised version [15] of the subtraction method [41]. 
The number of flavours was set to flve and the renormalisation (fiR) and factorisation 
(/ii?) scales were chosen to be fi^ = E^^-^t I^f = respectively. The strong cou- 
pling constant, a^, was calculated at two loops with A^g^ = 220 MeV, corresponding to 
as{Mz) = 0.1175. The calculations were performed using the MRST99 [8] parameterisa- 
tions of the proton PDFs. The jet algorithm described in Section |^ was also applied to 
the partons in the events generated by DISENT in order to compute the jet cross-section 
predictions. The results obtained with DISENT were cross-checked by using the program 
DISASTER++ [42]. The differences were always within 2% and typically smaller than 
1% [43]. 

Since the measurements refer to jets of hadrons, whereas the NLO QCD calculations 
refer to partons, the predictions were corrected to the hadron level using the MC models. 
The multiplicative correction factor (Chad) was deflned as the ratio of the cross section 
for jets of hadrons over that for jets of partons, estimated by using the MC programs 
described in Section |. In order to estimate the uncertainty in the simulation of the 
fragmentation process, events were also generated using the HERWIG 6.3 [44] program, 
where the hadronisation is simulated by using a cluster model [45]. The mean of the 
ratios obtained with ARIADNE, LEPTO-MEPS and HERWIG was taken as the value of 
Chad, since the three predictions were in good agreement. The value of Chad differs from 
unity by less than 10%, except in the backward region of the Breit frame where it differs 
by 20%. 

The NLO QCD predictions were also corrected for the Z-exchange contribution by using 
LEPTO. The multiplicative correction factor was deflned as the ratio of the cross section 
for jets of partons obtained with both photon and Z exchange over that obtained with 
photon exchange only. The correction is negligible for < 2000 GeV^ but reaches 17% 
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in the highest-Q^ region. 

Several sources of uncertainty in the theoretical predictions were considered: 

• the uncertainty on the NLO QCD calculations due to terms beyond NLO, estimated 
by varying fiji between and 2E:^-^^, was ~ ±5%; 

• the uncertainty on the NLO QCD calculations due to that on as{Mz) was estimated 
by repeating the calculations using two additional sets of proton PDFs, MRST99tt 
and MRST99II [8], determined assuming a>.(M^) = 0.1225 and 0.1125, respectively. 
The difference between the calculations using these sets and MRST99 was scaled by 
a factor of 60% to reflect the current uncertainty on the world average of Og [46]. The 
resulting uncertainty in the cross sections was ~ ±5%; 

• the variance of the hadronisation corrections as predicted by ARIADNE, LEPTO- 
MEPS and HERWIG was taken as the uncertainty in this correction, which was typi- 
cally less than 1%; 

• the uncertainty on the NLO QCD calculations due to the statistical and correlated 
systematic experimental uncertainties of each data set used in the determination of 
the proton PDFs was calculated, making use of the results of an analysis [10] that 
provided the covariance matrix of the fitted PDF parameters and the derivatives as 
a function of Bjorken x and Q^. The resulting uncertainty in the cross sections was 
typically 3%, reaching 5% in the high-E^^jg^. region. To estimate the uncertainties on 
the cross sections due to the theoretical uncertainties affecting the extraction of the 
proton PDFs, the calculation of all the differential cross sections was repeated using a 
number of different parameterisations obtained under different theoretical assumptions 
in the DGLAP fit [10]. This uncertainty in the cross sections was typically 3%. 

The total theoretical uncertainty was obtained by adding in quadrature the individual 
uncertainties listed above. 

6 Systematic uncertainties 

The following sources of systematic uncertainty were considered for the measured jet cross 
sections [43,47]: 

• the uncertainty in the absolute energy scale of the jets was estimated to be ±1% for 
-^Tjet > 10 GeV [48] and ±3% for lower -^rjet values. The resulting uncertainty was 
~ 5%; 

• the uncertainty in the absolute energy scale of the positron candidate was estimated 
to be ±1% [4]. The resulting uncertainty was less than 1%; 
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• the differences in the results obtained by using either ARIADNE or LEPTO-MEPS 
to correct the data for detector and QED effects were taken to represent systematic 
uncertainties. The uncertainty was typically smaller than 3%; 

• the analysis was repeated using an alternative technique [49] to select the scattered- 
positron candidate. The uncertainty was less than 2%; 

• the -E-Tjet cut was raised to 4 GeV. The uncertainty was smaller than 1%; 

• the cut in ?7j^^ used to suppress the contamination due to photons falsely identified 
as jets in the Breit frame was set to —3 and to —1.5. The uncertainty was typically 
~ 1%; 

• the uncertainty in the cross sections due to that in the simulation of the trigger and 
in the cuts used to select the data was typically less than 3%. 

In addition, there was an overall normalisation uncertainty of 1.6% from the luminosity 
determination, which was not considered in the cross-section calculation. 

The systematic uncertainties not associated with the absolute energy scale of the jets were 
added in quadrature to the statistical uncertainties and are shown on the figures as error 
bars. The uncertainty due to the absolute energy scale of the jets is shown separately as 
a shaded band in each figure, due to the large bin-to-bin correlation. 



7 Inclusive jet differential cross sections 

The differential inclusive jet cross sections were measured in the kinematic region > 
125 GeV^ and —0.7 < cos 7 < 0.5. These cross sections include every jet of hadrons in 
the event with E^^^^ > 8 GeV and — 2 < rj^^ < 1.8 and were corrected for detector and 
QED radiative effects. 

The measurements of the differential inclusive jet cross sections as functions of Q^, -E'^jet 
and ryj^t are presented in Figs. |l]— ^ and in Tables |l]— ^. The data points are plotted at 
the weighted mean in each bin of the corresponding variable as predicted by the NLO 
QCD calculation. The measured da/dQ"^ (dcr/dE^-^^) exhibits a steep fall-off over five 
(three) orders of magnitude in the (E^^^^) range considered. In the low-Q^ region 
(125 < Q"^ < 250 GeV^), the selected data sample covers 3 ■ 10~^ < x <2 ■ 10~^, whereas 
in the high-Q^ region (Q^ > 5000 GeV^), the range is 6 ■ 10^^ < x < 0.25. 

The measurements of the differential cross-section da/dE^-^^^ in different regions of 
are presented in Fig. ^ and in Tables § and |^. The -E^ jet dependence of the cross section 
becomes less steep as increases. 



8 



8 Comparison to NLO QCD calculations 



The NLO QCD predictions, corrected as described in Section ||, are displayed and com- 
pared to the measurements in Figs. jT]-^ It should be noted that the hadronisation cor- 
rection, shown in Figs. |I]c), |^c) and ^), was obtained with models (ARIADNE, LEPTO- 
MEPS and HERWIG) that implement higher-order contributions in an approximate way 
and, thus, their predictions do not constitute genuine fixed-order NLO QCD calculations. 
This procedure for applying hadronisation corrections to the NLO QCD calculations was 
verified by checking that the shapes of the calculated differential cross sections were well 
reproduced by the model predictions at the parton level. 

The ratios of the measured differential cross sections over the NLO QCD calculations are 
shown in Figs. |l|b), ||b), ^d) and |[ The calculations reasonably reproduce the measured 
differential cross sections, although they tend to be below the data. The agreement 
observed at high complements and extends an earlier comparison of the differential 
exclusive dijet cross sections at > 470 GeV^ [17]. For that measurement of the 
exclusive dijet cross sections, asymmetric cuts on the E^^^^t J^^^ were applied [17] 

to avoid infrared-sensitive regions where NLO QCD programs are not reliable [14]. This 
difficulty is not present in the calculations of inclusive jet cross sections and, as a result, 
the theoretical uncertainties are smaller than in the dijet case. Thus, measurements of 
inclusive jet cross sections allow more precise tests of the pQCD predictions than dijet 
production. 

At low and low E^^^ct^ calculations fall below the data by ~ 10%. The differences 
between the measurements and calculations are of the same size as the theoretical un- 
certainties. To study the scale dependence, NLO QCD calculations using ftR = fJ^F = Q, 
shown as the dashed line, are also compared to the data in Figs. |l|— ^; they provide a 
poorer description of the data than those using fiR = -E'l^jet- 

The overall description of the data by the NLO QCD calculations is sufficiently good to 
make a precise determination of a^. 

9 Measurement of ag 

The measured cross sections as a function of and -E^jct were used to determine 
as{Mz): 

• NLO QCD calculations of da/dA {A = were performed for the three 

MRST99 sets, central, TT and ii. The value of as{Mz) used in each partonic 
cross-section calculation was that associated with the corresponding set of PDFs; 
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• for each bin, i, in the variable A, the NLO QCD calculations, corrected for hadroni- 
sation effects, were used to parameterise the as{Mz) dependence of da/dA according 
to 

= Cl-as{Mz) + Ci,-as\Mz) ; (4) 

• the value of as{Mz) was then determined by a fit of Eq. (^) to the measured da/dA 
values for several regions of the variable A. 

This procedure correctly handles the complete as dependence of the NLO differential 
cross sections (the explicit dependence coming from the partonic cross sections and the 
implicit dependence coming from the PDFs) in the fit, while preserving the correlation 
between and the PDFs. 

The uncertainty on the extracted values of as{Mz) due to the experimental systematic 
uncertainties was evaluated by repeating the analysis above for each systematic check [43]. 
The overall normalisation uncertainty from the luminosity determination was also consid- 
ered. The largest contribution to the experimental uncertainty comes from the jet energy 
scale. 

The theoretical uncertainties, evaluated as described in Section ^, arising from terms 
beyond NLO, uncertainties in the proton PDFs and uncertainties in the hadronisation 
correction were considered. These resulted in uncertainties in as{Mz) of 3%, 1% and 0.2%, 
respectively. The total theoretical uncertainty was obtained by adding these uncertainties 
in quadrature. The results are presented in Tables ^ and 0. 

The best determination of as{Mz) was obtained by using the measured da/dQ'^ for > 
500 GeV^, for which both the theoretical and total uncertainties in as{Mz) are minimised. 
A good fit was obtained with = 2.1 for 4 data points. The fitted value is 

as{Mz) = 0.1212 ± 0.0017 (stat.)+S:E (syst.)^S:S (^h-) • 

As a cross check, the measurement was repeated using the five sets of proton PDFs of 
the CTEQ4 A-series [5]; the result is in good agreement with the above value. Two 
other determinations of as{Mz) were performed. The first made use of the measured 
da/dQ"^ for the entire Q"^ range studied, Q"^ > 125 GeV^, resulting in as{Mz) = 0.1244 ± 
0.0009 (stat.)^Q QQ^^ (syst.)^o'gQ4Q (th.). The second used the measured da/dE^^^^^ in 
the region where the hadronisation corrections are small, E^^^t > 14 GeV, resulting in 
as{Mz) = 0.1212 ± 0.0013 (stat.)+J;J|j^° (systO-S^SS (th-)- These results are consistent 
with the central value quoted above. 

The value of as{Mz) obtained is consistent with the current PDG value, as{Mz) = 
0.1181 ± 0.0020 [50] and recent determinations by the HI [21] and ZEUS [17, 19] Col- 
laborations. It is compatible with a recent determination from the measurement of 
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the inclusive jet cross section in pp collisions at y/s = 1800 GeV, as(Mz) = 0.1178 ± 
0.0001(stat.)^QgQg5(syst.)]'^QQQ75(th.) [51]. It is in agreement with, and has a precision 
comparable to, the most accurate value obtained in e~^e~ interactions [46]. 

The QCD prediction for the energy-scale dependence of the strong coupling constant has 
been tested by determining as from the measured differential cross sections at different 
scales. Since the NLO QCD calculations with fiR = E^^^et provide a better description of 
the data than those using fiji = Q, a. QCD fit to the measured da/dE^^^^^ was performed 
in each bin of -E'rjcf The principle of the fit is the same as outlined above, with the 
only difference being that the dependence of da/dE^-^^ in Eq. (H) was parameterised 
in terms of as{{E^^^^^)) rather than as{Mz), where {E^-^^^) is the mean value of -E'-fjet 
in each bin. The measured as{{E^ ^^^) values, with their experimental and theoretical 
systematic uncertainties estimated as for ^^(M^), are shown in Fig. ^ and in Table ||. The 
measurements are compared with the renormalisation group predictions obtained from the 
as{Mz) central value determined above and its associated uncertainty. The results are in 
good agreement with the predicted running of the strong coupling constant over a large 
range in E^-^^. 



10 Summary 

Measurements of the differential cross sections for inclusive jet production in neutral 
current deep inelastic e^p scattering at a centre-of-mass energy of 300 GeV have been 
presented. The cross sections refer to jets of hadrons identified with the longitudinally 
invariant kx cluster algorithm in the Breit frame. The cross sections are given in the 
kinematic region > 125 GeV^ and —0.7 < cos 7 < 0.5. 

NLO QCD calculations provide a good description of the measured differential cross 
sections for inclusive jet production at high Q^, > 500 GeV^, or high jet transverse 
energies, -E^ jet > 14 GeV. This observation complements and extends that of the exclusive 
dijet cross section to lower Q^. At low and low jet transverse energies, differences of 
~ 10% between data and calculations are observed, which are of the same size as the 
theoretical uncertainties. 

A QCD fit of the measured cross section as a function of for > 500 GeV^ yields 

as{Mz) = 0.1212 ±0.0017 (stat.)^™ (syst.)^SS (th-)- 

This value is in good agreement with the world average and is at least as precise as any 
other individual measurement. 
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g2 bin 
(GeV2) 



da/dQ^ 
(pb/GeV2) A 

stat 



^syst 



QED 
correction 


PAR to HAD 
correction 


0.950 


0.9283 ±0.0058 


0.947 


0.9463 ±0.0014 


0.959 


0.9542 ±0.0038 


0.955 


0.9579 ±0.0035 


0.963 


0.9623 ±0.0028 


0.918 


0.9727 ±0.0069 



125 - 250 
250 - 500 
500 - 1000 
1000 - 2000 
2000 - 5000 
5000 - 10^ 



1.107 
0.3714 
0.0919 
0.02068 
0.00325 
2.29 • 10^^ 



±0.018 
±0.0080 
±0.0029 
±0.00103 
±0.00024 
±0.40 • 10-^ 



+0.010 

-0.035 
+0.0038 

-0.0153 
+0.0008 

-0.0035 
+0.00055 
-0.00018 
+0.00021 
-0.00037 
+0.14 . 

-0.11 



+0.056 

-0.055 
+0.0156 
-0.0148 
+0.0031 
-0.0032 
+0.00047 
-0.00057 
+0.00004 
-0.00005 

+0.03 1 n-5 
-0.04 ■ 



Table 1: Inclusive jet cross-section da/dQ"^ for jets of hadrons in the Breit frame, 
selected with the longitudinally invariant cluster algorithm. The statistical, sys- 
tematic and jet- energy- scale uncertainties are shown separately. The multiplicative 
correction applied to correct for QED radiative effects and for hadronisation effects 
are shown in the last two columns. 



^Tjet bin 
(GeV) 


d(r/dE^.^, 
(pb/GeV) 


A Stat 


^syst 


Ajet-ES 


QED 
correction 


PAR to HAD 
correction 


8-10 


62.42 


±0.99 


+0.93 
-2.35 


+2.19 
-2.39 


0.955 


0.9170 ±0.0030 


10 - 14 


28.09 


±0.49 


+0.23 
-0.44 


+1.33 
-1.21 


0.951 


0.9488 ± 0.0033 


14 - 18 


10.66 


±0.29 


+0.05 
-0.39 


+0.49 
-0.53 


0.955 


0.9697 ±0.0039 


18 - 25 


3.16 


±0.12 


+0.04 
-0.15 


+0.17 
-0.14 


0.954 


0.9703 ± 0.0022 


25 - 35 


0.646 


±0.046 


+0.020 
-0.002 


+0.022 
-0.026 


0.944 


0.9698 ± 0.0026 


35 - 100 


0.0318 


±0.0043 


+0.0010 
-0.0023 


+0.0021 
-0.0014 


0.954 


0.9627 ±0.0082 



Table 2: Inclusive jet cross-section da/dE^-^^ for jets of hadrons in the Breit 
frame, selected with the longitudinally invariant cluster algorithm. For details, 
see the caption of Table 





bin 


da/dr]f^^ 
(pb) 


A Stat 


^syst 


^jct-ES 


QED 
correction 


PAR to HAD 
correction 


-2 - 


- -1 


5.27 


±0.36 


+0.09 
-0.16 


+0.25 
-0.21 


0.942 


0.798 ±0.016 


-1 - 


-0.25 


46.5 


±1.2 


+0.9 
-1.4 


+2.7 
-2.8 


0.947 


0.813 ±0.012 


-0.25 


- 0.25 


139.5 


±2.8 


+1.4 
-7.1 


+6.2 
-6.0 


0.953 


0.901 ±0.010 


0.25 


- 1 


157.7 


±2.7 


+1.0 
-3.9 


+6.1 
-6.1 


0.963 


0.9900 ±0.0040 


1 - 


1.8 


103.9 


±2.0 


+0.7 
-2.8 


+4.0 
-3.8 


0.957 


0.9982 ±0.0088 



Table 3: Inclusive jet cross-section da/drj^^ for jets of hadrons in the Breit 
frame, selected with the longitudinally invariant kx cluster algorithm. For details, 
see the caption of Table 



17 



£^T,jet bin 
(GeV) 


(pb/GeV) ^stat ^syst Ajct-ES 


QED 
correction 


PAR to HAD 
correction 


125 < < 250 GeV^ 


8-10 
10 - 14 
14 - 18 
18 - 25 
25 - 100 


32.90 ±0.71 ^^73 ^L38 

1 Q no _)_n QO +0.11 +0.73 
iO.UZ HZU.OZ _Q37 _Qgg 

3.75 ±0.16 t'of, 
0.895 ±0.059 l°:°f, 
0.0197 ±0.0027 t'oZl ^'oZl 


0.965 
0.963 
0.964 
0.963 
0.956 


0.9137 ±0.0029 
0.9380 ± 0.0075 
0.9496 ± 0.0069 
0.9394 ±0.0041 
0.9162 ±0.0067 


250 < < 500 GeV^ 


8-10 
10 - 14 
14 - 18 
18 - 25 
25 - 100 


1700 -Ln c;9 +0.27 +0.59 
11.00 HZU.OZ _Q54 _Qg5 

8.57 ±0.28 l',f, 

3.64 ±0.18 1^,1 
1 nn7 -i-n n79 +0.009 +0.068 

l.VVI ZSZV.VIZ, _o.092 -0.048 

0.0294 ±0.0036 l^S^^ ^^SJ^ 


0.949 
0.942 
0.953 
0.952 
0.937 


0.9205 ± 0.0084 
0.9573 ± 0.0043 
0.9748 ± 0.0064 
0.9685 ± 0.0042 
0.9539 ±0.0022 


500 <Q^ < 1000 GeV^ 


8-10 
10 - 14 
14 - 18 
18 - 25 
25 - 100 


7.83 ±0.36 tr,t t'oil 
3.77 ±0.18 trol t'oit 

1.87 ±0.13 1°:SI trol 

U./iO inU.UUZ „o.085 -0.029 

0271 ±0 0037 +°-°°" 

u.uz/i inu.uuo/ -0.0001 -0.0013 


0.938 
0.941 
0.949 
0.958 
0.951 


0.9205 ±0.0090 
0.9579 ± 0.0049 
0.9877 ±0.0041 
0.9888 ±0.0031 
0.9808 ± 0.0035 


1000 < < 2000 GeV^ 


8-10 
10 - 14 
14 - 18 
18 - 25 
25 - 100 


9 on -un 99 +0-02 +0.05 

Z.OU HZU.ZZ _QQg _QQy 

1.86 ±0.14 t'o'ol Tol 
1.006 ±0.099 t'oT, Tol^ 
0.287 ±0.037 t'ofi t'ofol 

0.0173 ±0.0030 1°::^! t'oZt 


0.934 
0.937 
0.945 
0.936 
0.943 


0.9170 ±0.0078 
0.9567 ±0.0041 
0.9856 ± 0.0049 
0.9976 ± 0.0006 
0.9946 ± 0.0038 



Table 4: Inclusive jet cross-section da/dE^^^^ in different regions of for jets 
of hadrons in the Breit frame, selected with the longitudinally invariant kx cluster 
algorithm. For details, see the caption of Table 
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^Tjet bin 
(GeV) 


(pb/GeV) '^stat ^syst ^jet-ES 


QED 
correction 


PAR to HAD 
correction 


2000 <Q^ < 5000 GeV^ 


8-10 
10 - 14 
14 - 18 
18 - 25 
25 - 100 


1.30 ±0.15 Tol Tol 
0.724 ±0.081 t^olt Tolt 
0.318 ±0.051 t'oZ 
0.209 ±0.035 1°:°°? l^o'ol 
0.0167 ±0.0032 t'oS ^VoZ 


0.934 
0.938 
0.941 
0.941 
0.950 


0.9143 ±0.0097 
0.9521 ±0.0058 
0.9869 ± 0.0076 
0.9955 ±0.0016 
1.0022 ±0.0039 


> 5000 GeV^ 


8-10 
10 - 14 
14 - 18 
18 - 25 
25 - 100 


258 ±0 073 +°-°29 +0-003 

U.ZOO ICU.UM _o.028 -0.009 

1 62 ±0 042 +°-°°^ +°-°°^ 

U.iUZ inU.U^Z _o.034 -0.003 
U.iiU ICU.UOZ _o.006 -0.005 

055 ±0 018 +°-°^^ +°-°°^ 
u.uoo inu.uio -0.000 -o.ooi 

OonS6 ±0 0014 +°°°°'^ +°-°°°^ 

U.UUOU ICU.UUi^ -0.0000 -0.0001 


0.998 
0.934 
0.937 
0.936 
0.927 


0.940 ±0.017 
0.958 ±0.012 
0.9777 ±0.0026 
0.9994 ± 0.0069 
1.00291 ±0.00090 



Table 5: Continuation of Table 0. For details, see the caption of Table 



region 
(GeV^) 


as{Mz) 


^stat 


^syst 




125 - 250 


0.1252 


±0.0013 


+0.0042 
-0.0048 


+0.0082 
-0.0062 


250 - 500 


0.1264 


±0.0017 


+0.0036 
-0.0047 


+0.0060 
-0.0036 


500 - 1000 


0.1203 


+0.0022 
-0.0023 


+0.0027 
-0.0038 


+0.0032 
-0.0016 


1000 - 2000 


0.1208 


±0.0032 


+0.0025 

-0.0021 


+0.0022 
-0.0022 


2000 - 5000 


0.1256 


+0.0047 
-0.0049 


+0.0045 
-0.0079 


+0.0040 
-0.0041 


5000 - 10^ 


0.1286 


+0.0146 
-0.0158 


+0.0055 
-0.0046 


+0.0045 
-0.0044 


> 125 


0.1244 


±0.0009 


+0.0034 
-0.0041 


+0.0057 
-0.0040 


> 500 


0.1212 


±0.0017 


+0.0023 
-0.0031 


+0.0028 
-0.0027 



Table 6: The as{Mz) values as determined from the QCD fit to the measured 
da/dQ"^, as well as those obtained by combining several regions in that distribution. 
The statistical, systematic and theoretical uncertainties are shown separately. 
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^Tjet region 
(GeV) 




^stat 


^syst 




8 - 


- 10 


0.1285 


±0.0015 


+0 0038 
-0.0052 


+0 0078 
-0.0046 


10 - 


- 14 


0.1238 


±0.0013 


+0.0037 
-0.0036 


+0.0060 
-0.0042 


14 - 


- 18 


0.1236 


±0.0017 


+0.0030 
-0.0039 


+0.0046 
-0.0035 


18 - 


- 25 


0.1188 


+0.0021 
-0.0022 


+0.0032 
-0.0037 


+0.0036 
-0.0026 


25 - 


- 35 


0.1157 


+0.0043 
-0.0044 


+0.0029 
-0.0027 


+0.0039 
-0.0029 


35 - 


100 


0.1422 


+0.0174 
-0.0178 


+0.0096 
-0.0112 


+0.0106 
-0.0088 


> 


8 


0.1241 


±0.0008 


+0.0034 
-0.0039 


+0.0055 
-0.0038 


> 14 


0.1212 


±0.0013 


+0.0030 
-0.0036 


+0.0041 
-0.0030 



Table 7: The as{Mz) values as determined from the QCD fit to the measured 
da/dE^ -^^j as well as those obtained by combining several regions in that distribu- 
tion. The statistical, systematic and theoretical uncertainties are shown separately. 



(^#.ei> 

(GeV) 




^stat 


^syst 




8.91 


0.2113 


+0.0042 
-0.0041 


+0.0109 
-0.0144 


+0.0228 
-0.0127 


11.65 


0.1851 


±0.0030 


+0.0087 
-0.0083 


+0.0141 
-0.0095 


15.70 


0.1721 


±0.0034 


+0.0059 
-0.0078 


+0.0094 
-0.0068 


20.69 


0.1538 


+0.0036 
-0.0037 


+0.0054 
-0.0063 


+0.0062 
-0.0044 


28.61 


0.1398 


+0.0064 
-0.0065 


+0.0044 
-0.0040 


+0.0057 
-0.0043 


41.98 


0.1660 


+0.0247 
-0.0240 


+0.0135 
-0.0153 


+0.0148 
-0.0121 



Table 8: The as{{E^-^^'j) values as determined from the QCD fit to the mea- 
sured da/dE:^-^^. The statistical, systematic and theoretical uncertainties are shown 
separately. 
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Figure 1: a) The differential cross-section da/dQ"^ for inclusive jet production 



with 



> 9> GeV and 



The inner error bars 



-2 < r/j^t < 1.8 (filled dots), 
represent the statistical uncertainty. The outer error bars show the statistical and 
systematic uncertainties, not associated with the uncertainty in the absolute energy 
scale of the jets, added in quadrature. The shaded band displays the uncertainty 
due to the absolute energy scale of the jets. The NLO QCD calculations, corrected 
for hadronisation effects and using the MRST99 parameterisations of the proton 
PDFs, are shown for two choices of the renormalisation scale, b) The ratio between 
the measured dajdQ'^ and the NLO QCD calculation; the hatched band displays 
the total theoretical uncertainty. The shaded band in c) shows the magnitude and 
the uncertainty of the parton-to-hadron correction used to correct the NLO QCD 
predictions. 
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Figure 2: a) The differential cross-section da/dE^-^^^ for inclusive jet production 
with -El^jet > 8 GeV and —2 < rj^^^ < 1.8 (filled dots). Other details are as described 
in the caption to Fig. |||. 
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Figure 3: a) The differential cross-section da/drj^^^ for inclusive jet production 
with -El^jot > 8 GeV and —2 < rj^^^ < 1.8 (filled dots). Other details are as described 
in the caption to Fig. 
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Figure 4: The differential cross-section da/dE^^^^ for inclusive jet production 
with -Ej^jct > 8 GeV and —2 < r/j^^ < 1.8 in different regions of (filled dots). 
Each cross section has been multiplied by the scale factor indicated in brackets to 
aid visibility. Other details are as described in the caption to Fig. 
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Figure 5: Ratios between the differential cross-sections da/dE^-^^ presented 
in Fig. ^ and NLO QCD calculations using the MRST99 parameterisations of the 
proton PDFs and 
caption to Fig. 



-^Tjct (fi^l^d dots). Other details are as described in the 
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Figure 6: The as{E^-^^^) values determined from the Q CD fit of the measured 
da/dE^^^^^ as a function of E^-^^. The inner error bars represent the statistical 
uncertainty of the data. The outer error bars show the statistical and systematic 
uncertainties added in quadrature. The dashed error bars display the theoretical 
uncertainties. The three curves indicate the renormalisation group predictions ob- 
tained from the as{Mz) central value determined in this analysis and its associated 
uncertainty. 
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